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Abstract 



Several models of baryo(lepto)- genesis that give rise to large inhomogeneities 
in the composition of the universe are presented. In particular, a variation of 
C*~) ■ primordial abundances of light elements by the factor 2-5 at large distances is 

predicted. A cosmological model of baryonic island is considered. Creation of 
domains with very large baryon and antibaryon number is described. Such domains 
mostly collapsed into primordial black holes that could be the dominant part 
of (non-standard) cosmological cold dark matter with a widely dispersed mass 
spectrum. A non-collapsed part of these bubbles might make clouds of matter or 
antimatter with an enriched abundances of metals. A mechanism for creation of 



^h. such exotic objects can be realized by mixed order phase transitions induced by 

r~| ! the infiaton field. 

> 

^ 1 Introduction 

It is commonly believed that the universe is the same everywhere (at least inside the 
present day horizon). This is the well known Copernicus or cosmological principle, or the 
principle of cosmic democracy. Observational astronomical data quite well agree with 
this principle and especially strong argument comes from the perfect angular smoothness 
of the cosmic microwave background radiation (CMBR). However, chemical inhomo- 
geneities may escape detection by CMBR either because they do not necessarily imply 
inhomogeneities in mass or energy densities or if their size is too small to be observable 
^lso: ITEP, Bol. Cheremushkinskaya 25, Moscow 113259, Russia. 
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by the contemporary CMBR telescopes. A detection of these objects by other astro- 
nomical methods may be inhibited by their large distance from the Earth. Thus at the 
present level of our knowledge, there seems to be plenty of room for "pieces" of universe 
with an exotic matter content. Of course the hypothesis of cosmologically large chemical 
inhomogeneities looks quite striking but it already existed in different forms for quite 
a long time (e.g. the idea of domains of cosmic antimatter). Several simple theoret- 
ical models predict such phenomena. The universe with non-homogeneous chemistry 
can be created if certain phase transitions occurred at inflationary stage leaving behind 
astronomically large bubbles with different physical conditions. Such phase transitions 
should have a rather unusual character. The effective mass of a scalar (Higgs-like) field 
might be negative but only during a finite (and rather small) period of time. This period 
should take place not too long before the end of inflation. These properties could be 
realized by a simple coupling of this Higgs-like field (order parameter) to the inflaton. 
The models with this property predict formation of chemically different domains in the 
early universe that might be astronomically large but not too large and sufficiently rare 
to escape existing observational bounds. 

Since this is the last school on cosmological phase transitions of this century and even 
of the Second Millennium it is tempting to speculate about cosmology in the coming 
XXI century. To verify or to disprove Copernicus principle seems to be a great challenge 
for the future (hopefully not too distant one). 

The content of the lecture is the following. First, we discuss the model of baryo(lepto)- 
genesis that could give rise to bubbles of chemically different phases in the universe. In 
sec. 3 we present a model of inhomogeneous lepto-genesis that predicts an existence 
of cosmologically large regions with strongly varying primordial abundances of light el- 
ements. The simplest version of the model predicts 2/3 of the sky with the so called 
normal abundances (25% of He 4 and about 3 • 10~ 5 of deuterium), 1/6 of the sky with 
much richer fraction of light elements (50% of He 4 and about 15 • 10~ 5 of D), and 1/6 of 
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the sky with a poor production of light elements (12% of He 4 and about 1.5 ■ 10~ 5 of D). 
In the next section an inhomogeneous baryogenesis leading to formation of baryonic and 
anti-baryonic cosmological islands is discussed. After that we consider a possible mecha- 
nisms of creation of cosmic antimatter in the universe dominated by baryons. Related to 
that is a model of formation of primordial black holes that could deliver a rather unusual 
form of cold dark matter with a very broad mass spectrum of constituent "particles". 

2 Baryo- and Lepto- Genesis 

After Sakharov [[Ij] formulated three general principles of baryogenesis: non-conservation 
of baryonic charge, breaking of particle-antiparticle symmetry, and deviation from ther- 
mal equilibrium, a variety of models were invented that predict, in particular, noticeable 
inhomogeneities in baryon or lepton number densities (isocurvature perturbations) which 
might lead to chemically inhomogeneous universe (for a review see ref. 0). We will de- 
scribe here one of the models that is especially interesting from this point of view, namely 
the Affleck and Dine scenario ||. The model is based on a supersymmetric theory that 
possesses the following generic properties. First, there exist scalar fields with non-zero 
baryonic and/or leptonic charges. We denote these scalar baryons or leptons generically 
as x- Such fields are obligatory in any supersymmetric model. For example they could 
be super-partners of quarks and leptons. Second, the potential of these fields x have the 
so called flat directions along which the energy does not rise. Such flat direction might 
exist, in particular, prior to spontaneous symmetry breaking before particle masses were 
generated by the Higgs mechanism. 

During inflation scalar fields (in contrast to fermions) may form a condensate by 
rising quantum fluctuations along flat directions ||, |5|, ^ |7j . A large baryonic charge can 
be stored in this condensate if baryon number is not conserved. The subsequent decay 
of this condensate in the course of processes that conserve baryonic charge would create 
a cosmological baryon (or lepton) asymmetry that might reach very large values, even 
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close to unity. Recall that the observed value is 

/3 = Nb/N^ « (3 - 5) • 1(T 10 (1) 

The model permits to create simultaneously a small baryon asymmetry in agreement 
with observations together with a large lepton asymmetry || [| . Moreover the asym- 
metries may be inhomogeneous and the scale of variation of lepton asymmetry may be 
much smaller then the scale of variation of baryon asymmetry || |2[ . 

To illustrate the main features of this scenario let us consider the toy model with the 
following potential term: 

U(x) = m 2 \ X \ 2 + Ailxl 4 + A 2 (x 4 + X* 4 ) (2) 

In the FRW cosmology, where the metric is given by ds 2 = dt 2 — a 2 (t)dr 2 , the equation 
of motion for \ has the form: 

(d 2 - d 2 /a 2 ) x + 3Hx + U\x) = (3) 

where H = a/a is the Hubble parameter. 

If A2 = Ai/2 the potential has flat directions in the limit m = and the field could 
easily evolve along them. It is convenient to introduce phase and absolute value of the 
field: 

X = re ie (4) 

and in terms of this quantities flat directions go along 1 + cos 48 = 0. 

The field x ma y possess leptonic or baryonic charge (or a mixture of both) and the 
corresponding vector current is given by 

3n = i (X*d^x ~ df,x*x) (5) 
According to the equation of motion (0) the current is not conserved if A2 7^ 0: 

d,f = \ X A - (x*) 4 l (6) 



In the homogeneous case the equation of motion takes the form: 

X + 3Hx + U'( X ) = (7) 

It is the Newtonian equation of motion for a point-like body in two dimensional space 
of the complex plane x m the potential U(x)- Here x is the acceleration, U' is the force, 
and 3Hx is the "liquid" friction term. Baryonic charge density is given by: 

B(t)=j = r 2 9 (8) 

and, in terms of the mechanical analogy, coincides with the angular momentum of the 
body moving in the potential U. If the potential is spherically symmetric, i.e. it depends 
only on r = \x\ the angular momentum (or baryonic charge) is conserved. The presence 
of valleys breaks the symmetry and accordingly the conservation of jo. 

During inflation, when/if m <C H, the field x is infrared unstable and its quantum 
fluctuations rise in accordance with the law @, ^|, 0: 

(X 2 ) = ^ - min {(- # 4 /m 2 ), (~ H 2 /V\)} (9) 

Hence the field x might acquire a large expectation value along flat directions and 
when inflation is over and a small slope of the potential created by m 2 |x| 2 becomes 
essential, the field returns to the origin releasing baryonic (or leptonic) charge by the 
decay into quarks or leptons. However, it is important to keep in mind that if the 
field evolved only in radial direction so that 9 = 0, the baryonic charge accumulated 
by quantum fluctuations would be zero, as is seen from eq. (§). Hence a motion in 
orthogonal direction is necessary for efficient baryogenesis. This motion may be produced 
by quantum fluctuations across the valley. So when the field x "lives" in the valley, far 
from the origin, it would oscillate between steep walls of the valley slowly approaching 
zero. When the field comes sufficiently close to the origin the quartic terms in the 
potential become negligible and the field evolution is governed by the quadratic term 
m 2 |x| 2 . In this region the field rotates around the origin with a decreasing amplitude. 
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The damping of oscillation both in the valleys and near origin is forced by the following 
two phenomena: by the Hubble friction that diminishes the amplitude according to 
\x\ ~ a~ 3 / 2 and by the particle production which is strongly model dependent and 
might be very much different for baryonic and leptonic x' s - Estimates of the particle 
production rate can be found in ref. 0. This effect may strongly damp motion along 
orthogonal direction in the valleys and significantly diminish the generation of the charge 
asymmetry. The decay of \ at this stage creates average zero baryonic charge because 
the angular velocity 9 frequently changes sign. On the other hand, the decay of the field 
on the final stage, when it rotates around the origin, proceeded with charge conservation 
and the produced quarks or leptons would have baryonic or leptonic charge equal to the 
initial magnitude of the angular momentum of \i a ^ the moment when \ arrived to the 
flat region and oscillations in angular direction changed into rotation. The magnitude 
of the latter is determined by the damping in the valley, while direction of rotation (i.e. 
the sign of baryon or lepton asymmetry) is determined stochastically by the direction 
of quantum fluctuations across the valley during inflation. In such a model C and CP 
non- conservation is not necessary. The latter is created by chaotic initial conditions. 
The size of domains with a certain sign of the charge asymmetry may be very large due 
to exponential expansion. However, a concrete realization of the model can meet serious 
problems and some fine-tuning of parameters or initial conditions may be necessary. 

The asymmetry generated by the decay of x a t the final stage could be very large if 
the life-time of % with respect to decay (or to the particle production rate) is sufficiently 
long. The oscillating massive field has the non-relativistic equation of state, P = (here 
P is the pressure density), so the energy density of such field decreases as 1/a 3 . It is 
one power of a slower than the decrease of the energy density of relativistic matter. The 
slower red-shift may create the situation when the matter in the universe is completely 
dominated by baryons (or leptons) with a negligible amount of antiparticles. 

If C (CP) violation is explicit then the picture may be somewhat different. In our toy 
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model we can realize this case assuming that coupling constant A and mass are complex, 
so that the potential takes the form: 

U( X ) = mV + K) V) 2 + Ax 4 + a*( x *) 4 (10) 

Both terms are necessary, otherwise the CP-odd phase can be rotated away. In this 
potential the field x{t) rnay acquire a non- vanishing angular momentum after it traveled 
down along the bottom of the valley without any orthogonal motion and closer to the 
origin it changes "rails" from the lamb da- valley to the mass- valley that has a different 
direction. The transition from the bottom of one valley to another creates a transverse 
motion or, in other words, a nonzero baryonic or leptonic charge. 

3 Large and Inhomogeneous Lepton Asymmetry 

Models based on the picture described in the previous section allow to obtain a small 
baryonic and a large leptonic asymmetry if e.g. the decay of the baryonic x-held is 
faster than the decay of the leptonic one. In this case baryonic angular fluctuations in 
the valleys would be stronger suppressed than leptonic ones, so the initial value of the 
baryonic charge at the beginning of rotation would be smaller than the initial value of 
the leptonic charge. Moreover, the dilution of the produced baryons by expansion could 
be stronger than the dilution of later created leptons. 

The characteristic wave length of the variation of the baryon (Is) or lepton (II) 
charge asymmetry was estimated in ref. || and is approximately equal to the smaller of 
H^ 1 exp(l/vX) and Hf 1 exp[(ifj/m) 2 ]. In particular, it is possible to create the universe 
with a small and homogeneous baryon asymmetry and large (of the order of unity) 
lepton asymmetry, which in addition may be inhomogeneous at the scales well inside 



the present day horizon || |TD, [11] . Such domains of leptonic charge would create very 
strong density inhomogeneities and, in turn, too large angular fluctuations of CMBR. 
However the flavor symmetry between different lepton families, e <-> «-> r, permits to 
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avoid this problem because in this case one would expect domains that are symmetric 
with respect to permutation of different leptonic charges and in first approximation the 



density contrast is vanishing [10, 11 



In the simplest version of this scenario the universe would be filled with equally prob- 
able domains with (L e , L^, L T ) = (±1, 0, 0) and with all possible permutations between 
leptonic charges. Evidently there are no walls between such domains, as follows from 
the way of their creation, and so the cosmological problem of heavy domain walls |12] is 
avoided. Since the abundances of light elements produced at big bang nucleosynthesis 
are especially sensitive to charge asymmetry of electronic neutrinos, there would be two 
types of anomalous domains with low and high abundances and 4 types of the normal 
ones with the same normal abundances, as is written in the introduction. 

As discussed in the papers [[K], [ll]] the characteristic distances between the domains 
should be larger than a few hundred Mpc to agree with the observed small angular fluc- 
tuations of the CMBR temperature, 5T/T. Another possibility of small size domains, 
that could create a large 5T/T on very small scales, which are allowed by observations, 
was not explored but possibly, though not surely, this option is excluded by direct ob- 
servation of the light element abundances. Still it remains surprising that twice larger 
mass fraction of He A at not so large distances is not yet observationally excluded. An- 
other signature of such inhomogeneities in primordial He A is a strong variation of the 
exponential slope in 5T/T at small angles (at / > 10 3 ) at different patches on the sky. 
Possibly it is the most promising way to obtain an upper limit on inhomogeneities of 
primordial He 4 . 

4 Inflaton Induced Phase Transitions and Baryonic 
Islands 

If the field that is responsible for baryo(lepto)-genesis is coupled to inflaton then a 
very interesting pattern of charge asymmetries may be created. This field <fi could be 
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e.g. the Affleck-Dine one described in the previous section or the complex scalar field 
whose condensate creates spontaneous breaking of charge symmetry ||13|| . The coupling 
to inflaton, $, is assumed to be of the general renormalizable form: 

C mt = A|0| 2 $ 2 + s|0| 2 $ = A($ - ^) 2 \<f>\ 2 - A$ 2 |0| 2 (11) 

where $1 = g/2X is a constant. We assumed that the interaction conserves C (CP), 
though it is not necessary. 

Generically the effective mass of the field <fi may contain the following contributions: 

(mi ff f =ml + {R + (3T 2 + A($ - $x) 2 (12) 

where R is the curvature scalar and the correspoding term comes from a possible non- 
minimal coupling to gravity. The third term comes from the interaction with thermal 
bath, and the last one comes from the interaction with inflaton (pTTJ) . Note that at 
inflationary stage R = 12H 2 , while at radiation dominated stage R = 0. With such a 
form of mass term one can easily visualize the following picture. At an early inflationary 
stage the effective mass squared of <fi is large and positive, so that (f) sits near the origin 
= 0. When the inflaton evolves down and approaches the value $ ~ $i the mass 
squared becomes negative and remains such for some period of time determined by the 
inflaton evolution. When the amplitude of the inflaton tends to zero and inflation ends, 
m 2 jj may become positive again. During a relatively short period of negative m 2 the 
road is open for <fi to travel away from the origin and its subsequent destiny depends 
upon its self-interaction potential. The field may either return to the origin, or stay 
for a while in a false vacuum state with ((f)) ^ 0, or evolve far away along one or other 
flat directions of the potential. Correspondingly the generation of baryon (or lepton) 
asymmetry would be quite different. Some striking possibilities are discussed in the 
papers |L4], [T5|, [16j or in the review 0. 

In particular, our universe may be a huge baryonic island floating in the sea of dark 
matter fH] , This can be achieved if the field <ft is responsible for the spontaneous 
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breaking of C and CP symmetries so that cosmological charge asymmetry might 



~UJ (13) 



be generated only in the regions where (<p) 7^ 0. The size of this island may be of the 
order of the present day horizon or much larger. In the first case a small observed dipole 
asymmetry of CMBR forces us to live very close to the center of the island and it makes 
the model rather unnatural. Its interesting property is that the angular distribution of 
the CMBR temperature in this model would have an intrinsic dipole moment: 

'ST\ ( d 

dipole 

where Ri is the radius of the island and d is our displacement from the Center of the 
world. In this case the dipole contribution into angular fluctuations of CMBR is partly 
given by our peculiar motion and partly by the displacement from the Center and one 
does not need the very large value of the velocity, about 600 km/sec, whose origin is 
mysterious, to explain the observed dipole. There should also be an intrinsic quadrupole 
with a rather large value Q ~ (d/Ri) 2 , which is close to the observed one. 

The size of the island may be of super-horizon scale but nevertheless the angular 
fluctuations imprinted on CMBR by the existence of the boundary could be observ- 



able [ 21 1 . Of course their magnitude would be much smaller than in the sub-horizon 
case and we may live far from the Center without breaking the observational limits. 
Still the large magnitude of the dipole asymmetry of CMBR may be explained by the 
insular hypothesis. 

Models of baryogenesis with spontaneous breaking of charge symmetry predict for- 
mation of cosmological matter-antimatter domains, as was noticed in ref. [|17|]. The size 
of these domains may be cosmologically large if after their formation the universe passed 
through a period of exponential expansion (inflation) ]18| . A review of the earlier ideas 
on the subject can be found in the paper |1| or in the recent one [2(J, which contains a 



discussion and the list of references to later development. In the model described above 
C(CP) is broken spontaneously but only in finite bubbles, while outside the bubbles par- 
ticles and antiparticles are completely symmetric and cosmological charge asymmetry 
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in these regions is vanishing. As in any model with spontaneous breaking of charge in- 
variance the universe as a whole should be charge symmetric, so there must be an equal 
number of bubbles and anti-bubbles. We assumed above that the size of the bubbles 
(or at least of our bubble) is very large, close to the horizon or much larger than it. 
In such a model other bubbles most probably should be far beyond the horizon and be 
unobservable. 

In principle another scenario is possible when the bubbles are relatively small so that 
there are many of them inside the horizon and the universe reminds a piece of Swiss 
cheese where bubbles of matter and antimatter are separated by voids filled with dark 
matter. There is no direct contact between them and thus no pp-annihilation. This weak- 



ens the limits on antimatter domains found in refs. p2| , P3d but the observed distribution 
of baryonic (luminous) matter seems to disagree with such a picture. Moreover, angu- 
lar variation of the CMBR temperature induced by such isocurvature perturbations are 



likely to be unacceptably large. A related analysis was performed in references |24], |25 
but the case considered there rather differs from this model and some more work is 
desirable. 



5 Large Baryon Asymmetry, Primordial Black Holes, 
and Dark Matter 

In the model of the previous section the magnitude of baryon asymmetry inside the 
islands was normal (]l|), but the sign could be different. The sea outside the islands 
was charge symmetric. Here we will consider another exotic scenario of baryogenesis 
according to which the average baryon asymmetry is normal and positive but in its 
background there could be relatively small (but still astronomically large) bubbles with 
a very large (even close to unity) asymmetry. Such a picture can be realized with the 
help of the scalar baryon discussed in sec. |2| if the flat directions of its potential are 
"locked" by a positive value of m^ff during almost all inflationary stage and the window 
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is only open when the inflaton field $ is close to $1, see eq. (0). In such a scenario the 
Affleck-Dine field \ may reach large values and accumulate a big baryonic charge only 
in relatively small bubbles, where the field succeeded to reach a certain critical value. 
This process is stochastic and its probability may be sufficiently small so only a minor 
fraction of the universe would be occupied by this baryon rich bubbles (see ref. Jl6| for 
the details). If the potential of is C(CP)-symmetric then the probability of baryon 
rich and antibaryon rich bubbles would be the same. So according to this scenario 
the cosmological background is baryo-asymmetric and each bubble is (strongly) baryo- 
asymmetric but on the average there is equal amount of baryonic and anti-baryonic 
bubbles. 

The subsequent destiny of these baryon rich bubbles depends upon their size and the 
value of the baryon asymmetry. Predominantly they would end as primordial black holes 
formed at the stage when baryons became nonrelativistic. According to calculations of 



ref. [16 1 the mass spectrum of the black holes in a simple version of the model is given 



by the log-normal distribution: 



dN n ( l 2 M ' 

Cexp -7 In 2 — (14) 



dM r V Mi. 

where C, 7, and M\ are constants expressed through some unknown parameters of the 
underlying theory. We cannot make any reliable prediction about their value. They may 
easily vary by several orders of magnitude. So in what follows we will treat them as free 
parameters to be constrained from observations. 

According to this model not only the bubble sizes are stochastically distributed but 
also the magnitude of baryon asymmetry inside them. Depending upon the values of 
these two parameters the bubbles would either form black holes and whether they are 
baryonic or anti-baryonic could not be distinguished, because black holes do not have 
baryonic hairs. Smaller size objects or those with a smaller value of the asymmetry 
do not collapse and on the tail of the distribution we may expect to observe clouds of 
antimatter in close vicinity, anti-stars, or even possibly small anti-galaxies. 
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Another interesting feature of this baryon (or antibaryon) rich regions is that their 
primordial chemistry should be very much different from the standard one. The out- 
come of light elements produced during "the first three minutes" strongly depends upon 
the baryon-to-photon ratio, Nb/N 1 . Normally it is very small, eq. ([!]), and the nucle- 
osynthesis effectively stops at He 4 . The amount of the next stable light element Li is 
roughly 9 orders of magnitude smaller. In the case of a larger ratio Nb/N 1 the amount 
of heavier elements could be much larger and the baryon rich bubbles might have much 
more evolved chemistry which is more typical for the system contaminated by stellar 
processes. 

The described here mechanism easily explains an early production of quasars, while 
in the standard theory of structure formation this creates a serious problem. The model 
also explains an evolved chemistry around quasars. This could be the matter from 
baryon rich region that did not completely disappeared in the black hole (quasar) and 
which was enriched by heavier elements from the "very beginning", i.e. from big bang 
nucleosynthesis. In the standard approach it is assumed that heavier elements observed 
in the vicinity of quasars were produced as a result of stellar processes in the stars of first 
generation. But this is very difficult to accomodate to the standard theory of structure 
formation according to which structure are formed much later than the observed distant 
quasars. 

An interesting implication of the model of primordial black hole formation with the 
mass spectrum fll4|) is that they might make a new form of cold dark matter. In contrast 
to the usual type of dark matter, when all dark matter particles have equal masses, this 
model predicts a mass dispersed CDM where particles of different masses play different 
role. The parameters describing the spectrum can be fixed from the observations. In 
particular, one could request that there is one super-heavy black hole (M > 10 6 M Q ) per 
large galaxy. Here M & is the solar mass. These super-heavy black holes can serve as seeds 
for galaxy formation. The bulk of dark matter may be presented by black holes with 
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the masses near 10 — 100 M , while those with the mass about 0.5 M Q may explain the 



MACHO 's observations [pq , |27l- In addition the model gives an explanation of formation 
of heavy (M > 100M Q ) black holes observed in galaxies. A possible simple fit: 7 = 0.14 
and Mi = O.5M , gives for the relative contribution to the universe mass density in 
different mass intervals the following results: F(M < 0.1) = 0.3%, F(0.1 - 2) = 12%, 
F(M < 10) = 38%, F(M = 2 - 30) = 48.5%, F(M = 30 - 1000) = 37.5%, F(M = 
10 3 - 10 6 ) = 1.6%, F(M = 10 6 - 10 8 ) = 3.6 • 10~ 9 , where M is expressed in solar mass 
units. Of course a more accurate analysis with the account of existing astronomical 
bounds on the cosmic density of black holes in different mass intervals is necessary to 



make more definitive conclusion and numerical simulation of structure formation 
in CDM model with the mass spectrum of dark matter particles given by eq. fli~4]) is 
desirable. 
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